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ABSTRACT 



Accurate ocean data is essential for successful fleet operation. The N-ROSS 
Satellite, which is being developed for this mission, will carry a Low Frequency 
Microwave Radiometer ( LFMR ). The LFMR consists of large flexible reflector and 
boom and spins at 15 r.p.m. The effects of the flexibility of the boom, the spin-up 
procedure and the structural damping on the pointing error of the LFMR are 
investigated by performing the dynamic simulation using the Dynamic Simulation 
Language. Two cases of boom material. Aluminum Alloy and the Graphite.epoxv 
composite material, are analyzed and the results are compared. The simulation and 
analysis results are presented in graphical forms. 
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I. INTRODUCTION 



A. BACKGROUND 

Accurate ocean weather prediction is essential for successful fleet operations. The 
NANY needs superior data collection capability to obtain the data density and 
reliability necessary to produce consistently accurate forecasts and oceanographic data. 
However present prediction models are limited by the quality and quantity of input 
data which come mainly from ships. The obvious approach to satisfy these purpose 
can be derived from satellite observations. Therefore, the NAVY planned the 
construction of the Navy Remote Ocean Sensing System ( N - ROSS ). [Ref. 1] This 
system consists of satellites which scan the earth surface and provide the fleet with 
timely worldwide knowledge of ocean data such as seasurface wind speed, wind 
direction, seasurface temperature, ice edge detection, ocean wave height and ocean 
photograpy. [Ref. 2] To satisfy the mission requirements the N - ROSS Satellite will 
carry several sensors. Among these sensors the Low Frequency Microwave Radiometer 
( LFMR ) is the most important and the most interesting from the dynamics of the 
spacecraft point of view. The function of the LFMR is scanning the earth surface and 
measure the seasurface temperature. To increase the scanning area the deployable 
reflector spins at 15 r.p.m. The sizes of this LFMR reflector and boom are relatively 
large compared to the N - ROSS Satellite itself. So the weight of this boom should be 
light, which makes the boom flexible. By this reason, there exist certain extent of 
deflection at the tip of the LFMR boom and this boom vibrates when this boom is 
spinning. Deflection and vibration due to this elastic deformation induce pointing 
error of the reflector in elevation and azimuth angle. However there is strict pointing 
error requirement of the LFMR. [Ref. 2] Therefore analysis of the flexible LFMR 
boom which supports the sensor payload is imperative for this research. 

B. STATEMENT OF PROBLEM 

The traditional approach to dynamics of a boom system is based on the 
assumption that the systems are composed of rigid bodies. Until recently, only a rigid 
body motion was assumed for the analysis. However, the flexible system includes a 
small elastic deformation as well as a large motion. These small elastic deformations 
include bending, twisting and axial extension. Development of a dynamic model 
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including flexibility demands more accuracy for the system responses. Without 
considering these small motion of a boom, we cannot expect a certain accuracy to 
maintain a spacecraft attitude and pointing control. [Ref. 3] Recently, efforts have 
been made to control maneuvers of mechanical systems which can not be adequately 
modeled using a rigid body assumption for all or some of the system components, 
especially in the fields of satellites, [Ref. 4: pp. 257-264] 

Therefore, the development of a good dynamic model of a flexible system, an 
efficient dynamic equations formulation method and a good dynamic simulation schem 
are essential for the analysis of and identification of potential problems in the flexible 
LFMR system. 

C. THESIS OUTLINE 

In Chapter II, the development of an analytic model for the Lower Frequency 
Microwave Radiometer ( LFMR ) reflector boom in 3-dimensional motion is described. 
The large motion due to rotation is described by an equivalent rigid boom motion and 
elastic deflection of a flexible boom relative to the equivalent rigid boom motion is 
expressed using the mode superposition technique. The dynamical equations for this 
model are formulated using the Lagrange's method. 

In Chapter III, The computer implementations for the solution of the obtained 
equations are explained. For the modal analysis of the system. NAsa STRuctural 
ANalysis ( NASTRAN ) computer program was used and Dynamic Simulation 
Language ( DSL ) was applied to solve the simultaneous, nonlinear, ordinary 
differential equations. The L INPACK subroutines DGEFA and DGESL are also used 
in the dynamic simulation. 

In Chapter IV, simulation results are presented to investigate the deflection and 
pointing error of the LFMR. Comparisons are made by changing the torque input 
condition. The problems considered are 1) the effects of spin-up procedure on the 
pointing error of the LFMR reflector; 2) the effect of damping on the settling time of 
pointing error, 3) the equilibrium configuration of LFMR booms due to constant 
rotating speeds For the comparison purpose, two kinds of material, aluminum alloy 
and composite material, were assumed as the LFMR boom materials. 

In chapter V, Conclusions are made from the research and some 
recommendations for future work in the area of the dynamic analysis of the flexible 
LFMR reflector boom system are given. 



12 



II. FORMULATIONS OF DYNAMIC EQUATIONS 



A. INTRODUCTION 

In this Chapter, the dynamical equations of motion for the LFMR system which 
rotates in three dimensional space is developed. The 2-dimensional planar motion of 
the same boom was also studies and are presented in Appendix A. A simple dynamic 
model of the LFMR system is developed for this analysis. The Lagrangian approach 
and the mode superposition technique are used for the formulation of the dynamic 
equations of the flexible LFMR system. 

B. DESCRIPTIONS OF THE MODEL 

The LFMR shown in Fig. 2.1 consists of four structures: a reflector, an upper 
reflector boom, a lower reflector boom and an electronics box. The reflector is attached 
to the top of the upper reflector boom. The upper reflector boom and the lower 
reflector boom is connected by a boom hinge. The electronics box is attached to the 
bottom end of the lower boom. 

For our analysis, the deployable reflector is modeled as a concentrated mass at 
the tip of the upper reflector. We assume the boom hinge which connects the two 
booms is stiff and firm and there is no relative motion between the booms after the 
deployment of the LFMR boom. Therefore we consider the whole system ( reflector, 
booms, boom hinge ) as a one body system. 

The LFMR system is connected to the Main Bus of the X - ROSS Satellite by a 
Spacecraft Boom. The attitude of the X - ROSS Satellite is controlled by a Attitude 
Determination And Control System ( ADACS ) very accurately. [Ref. 5] Therefore, it 
is assumed that the spin axis and the base is remain fixed in the reference frame fixed 
to the X - ROSS Satellite. The X - ROSS Spacecraft moves on a circular orbit with 
the spin axis always pointing the earth center. Hence, the gravitational force is in 
equilibrium with the centrifugal force in the orbit plane: the LFMR system is in zero-g 
environment. Therefore, in the dynamic model of the present studies, the reference 
frame fixed to the X - ROSS Spacecraft is assumed the Xewtonian ( inertial ) reference 
frame and the LFMR system is in zero-g environment. 

From the above assumptions the dynamic model of the LFMR system is defined 
as shown in Figure 2.2. The global coordinates X. Y, Z is fixed in the inertial reference 
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frame and a moving coordinate x, y, z ( local coordinates ), which is a body fixed 
coordinate system, is defined as shown in Figure 2.2. The body fixed coordinate 
system ( local coordinates ) is attached to the base O. The local z-axes and the global 
Z-axes are the common axis of rotation. The local x, y-axes and the global X, Y-axes 
are in the same plane with angle difference 0. 

C. LAGRANGE'S EQUATION 

For any system there must be same numbers of independent coordinates as the degrees 
of the freedom of the system to completely describe the motion of the system. The 
choice of coordinates is important in dynamic analysis. Such independent coordinates 
are called generalized coordinates and are denoted by the letter q f . For a system with a 

set of n independent generalized coordinates q r ( r = 1, 2. 3 n ), Lagrange's 

equations are expressed as [Ref. &7] 



d 

dt 




<3T <3U 

-r 

dq k oq k 



(eqn 2.1) 



(k = 1,2.3, , n ) 

where T is the kinetic energy, U is the potential energy and is the generalized forces 
which is defined as follows 



_ 6R. 

Q. = V F. (eqn 2.2) 

- 1 n 

The dot over a variable means derivative with respect to time. F is the force acting on 
particle j and R : is the instantaneous position of particle j and may be expressed in 
terms of generalized coordinates 

Rj = R. ( q v q 2 , q ? , q n ) (eqn 2.3) 

and 5Rj is the virtual displacement of the particle. 

To describe the motion of the LFMR boom system which composed of a large 
slow motion due to rotation and a small fast motion due to elastic vibration, two kinds 
of generalized coordinates are defined. One is 9 for rotation of boom and others are q h 
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PER REFLECTOR 




Figure 2.1 N-ROSS Baseline Configuration. 
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M 




Figure 2.2 Modeling of LFMR boom system. 

( h = 1, 2. 3 n ) for h-th mode generalized displacement, where n is number of 

modes. 

Now Lagrange's equation 2.1 is rewritten as 



d 

dt 



dT 

W ] " 



dT 

4 - 

(50 



d\J 

<50 



O 



(eqn 2.4) 
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and 



d 



dt 




8T 8\J 

+ = Qu 

h 



(eqn 2.5) 



( h = 1, 2, 3, , n) 

The external force acting on the LFMR system is assumed the torque t by a 
torque motor at the root of the boom. The contribution of this torque to the 
generalized forces is Q r = t and this torque does not contribute to Qj . Damping 
forces are assumed equal to the modal damping value in this analysis. Since the modal 
damping values can be measured easily. Therefore the contribution of damping forces 
to the generalized forces can be obtained using a dissipation function. 

D = — l — X 2 ^ to i Mj qf (t) (eqn 2.6) 

5D 

= “ 2 ^h 03 h M h % 

where 

: modal damping ratio of h-th mode 

to : the natural frequency of h-th mode 

: the modal mass of h-th mode 

Then equation 2.4 and 2.5 finally written as 



with 



d 

dt 



31 

~3Q 



3 T 8 U 

<30 <30 



and 





<3t au 

d % 




(0 h M h % 



(eqn 2.7) 



(eqn 2.8) 
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( h = 1,2,3, 



n ) 

D. POSITION AND VELOCITY 

During the rotational motion of the LFMR system, the boom deforms. The 
deformed positions of a generic point in the system can be expressed the vector sum of 
a vector Rq(x) from the origin to the undeformed position of the point and a vector 
W(x.t) which is deflection, as shown in Figure 2.3. The notations show in the Figure 
2.3 represent the following parameters: 



M : 


tip mass 


n V : 


mass of electronics box 




length of lower boom 


t 2 : 


length of upper boom 


P*: 


angle between two links 


Z : 


applied torque 


0(0 : 


angular displacement 


0(t) : 


angular velocity 


R 0 (x) : 


position vector of the point on the 


R(x.t) : 


position vector of the point on the 


W(x,t) : 


deformation vector of boom 


i : 


unit vector of local x-direction 


j : 


unit vector of local v-direction 


k : 


unit vector of local z-direction 


j o : 


unit vector of global X-direction 


: 


unit vector of global Y-direction 


kg : 


unit vector of global Z-direction 



Then the position vector R(x.t) is expressed as 



R(x.t) = Rq(x) -r W(x.t) (eqn 2.9) 

The undeformed position Rg(x)is represented by its components, 

R 0 (x) = R„(x) + R z (x) (eqn 2.10) 

= R x (x) i + R z (x) k 
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Figure 2.3 Parameters of the boom system. 

The elastic deflection W(x,t) is expressed as the modal sum as follows: 

W(x.t) = V [ 9j X (x) i + (pTlx) j + (p/(x) k ] qj(t) (eqn 2.11) 

i 

where 
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(pj X (x) i is i-'.h mode shape function in extension 
<p ; y (x) j is i-ih mode shape function in translation 

From the equation 2.10 and equation 2.11, equation 2.9 can be rewritten in the form of 

R(x,t) = R x (x) i + R z (x) k (eqn 2.12) 

+ 1 [ <P; X (x) i + <P-, y ( x ) j + <p; z (x) k ] q.(t) 
i 

= [ R x (x) + (p x (x) q ; (t) ] i + [ £ <p y (x) q.(t) ] j 

i 

+ [ R z (x) + I q>i Z (x) q.(t) ] k 
i 

The velocity of the point in the Newtonian reference frame is -obtained by taking the 
time derivative of equation 2.12 and applying the relation between the time derivative 
of unit vector i and j, 

i = 0 k x i = 0 j 

j = 0 k x j = — 0 i 
k = 0 k x k = o 

the velocity is 



R(x.t) = [ - 0 y <p; y (x) q ; (t) + X ( Pf X ( x) ^(0 1 ‘ ( ec l n 2 - 13 ) 

i i 

+ [ 0 R x (x)+ 0 1 <p-, x (x) q ; (t) + 1 <Pi y ( x > ^(0 1 i 
i i 

+ [ X <Pi Z ( x ) q/t) ] k 
i 



E. KINETIC ENERGY 

The kinetic energy of the system can be expressed by the summation of three 
different kinetic energies. One is the kinetic energy of the boom itself, another is the 
kinetic energy of the tip mass and the other is the kinetic energy of the R.F electronic 
box attached to the origin of the boom. 
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The kinetic energy of the boom itself is 



T 



bm 



i f e. 

— f R(x,t) • R(x,t) dm 
— J* R(x,t) • R(x,t) pdx 



where 



dm : 



dx : 

P : 

e : 



R(x,t) : 



differential mass of the boom 

differential length of the boom 

mass per unit length 

total length of the boom 

absolute velocity of the points on the boom 



The kinetic energy of the tip mass is 



T = 

tm 



M R(C,t) • R(C.t) 



where 

M : tip mass 

R(£,£) : the absolute velocity at the tip position 

The kinetic energy of the R.F electronic box is 



T rf = 



1 

T'r 



e 2 



zz 



where 



I- : mass moment of inertia for the electronic box 

• l zz 

0 : time derivative of aneular displacement 



Thus the total kinetic energy of the system is 



(eqn 2.14) 



(eqn 2.15) 



(eqn 2.16) 
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(eqn 2. 17) 



T " T bin + T isn + T rf 



i f e . 

— J R(x,t) • R(x,t) dm 



M R(C.t) • R(C.t) 



+ — I 0 2 
2 r zz 



F. POTENTIAL ENERGY 

The potential energy U of the flexible boom system can be composed of the 
gravitational potential energy L* CT due to rotation of the system and the strain energy 
(or elastic energy) L’ s of the boom due to deformation. But in our model analysis, we 
exclude gravitational acceleration and only consider the strain energy. The potential 
energy was determined by the work done by the static weight in the deflection. This 
work is, of course . stored in the flexible flexible boom system as strain energy. 

In this thesis, we apply mide summation method to expand the deflection in terms of 
the normal modes of the system. The deflection of a boom without any external forces, 
satisfies the lollowing equation of motion. [Ref. 7] 

[ EIW"(x,t) ]" + p\V(x,t) = 0 (eqn 2.18) 

where 

p : mass per unit length 

El : flexural rigidity 

and ' represent the derivative with respect to x. The normal modes (p;(x) of the boom 
must satisfy the equation 

[ EI<pf(x) ]" - co-“ p <pj(x) = 0 (eqn 2.19) 
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and its boundary conditions. The normal modes (pj(x) also satisfy the orthogonality 
relation 



t 

j^(p-(x) cpj(x) dm' = 0 (for i*i) (eqn 2.20) 

= M i (for j = i ) 

where VI • is the generalized modal mass of the i-ih mode. 

As expressed in Appendix A, the deflection of the boom in the equation A. 4 for 
the general form is 

W(x,t) = 2>i(x) qi (t) (eqn 2.21) 

i 

and the generalized coordinate q'j(t)'can be determined by applying Lagrange's 
equation after setting up the kinetic and potential energies. 

Now the potential energy can be expressed as 



U = 



1 r C 2 

— f EIW" 2 (x,t) dx 

2 J 0 



y-£ I % 9j Jl EIcpj'Xx) (Pj"(x) dx 

i j 



(eqn 2.22) 



Vlultiplying the both sides of equation 2.19 by <Pj(x) and after integration for the whole 
boom, equation 2.19 becomes 



J 



6 

( ^(pj(x) [ El <pf(x) ] dx 



®r J 0 ( pi( x ) ^j( x ) dm 



(eqn 2.23) 



After integration by parts and using the boundary conditions, equation 2.23 becomes 



f 



e 

o EI [ <pf(x) cpj"(x) J dx 



(eqn 2.24) 
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= 0) i J 0 <Pi(x) <pj(x) dm 



From the orthogonaiity condition, the equation 2.20 and from equation 2.24, the strain 
energy in equation 2.22 becomes 

u = -i- Z C 0 £ 2 M: qi (tr (eqn 2.25) 

z i 

G. DERIVATIONS OF EQUATIONS OF MOTION 

From the equation 2.13. the dot product of R(x,t) is 

R(x.t)*R(x,t) = [ - 6 V (p.- v (x) qi (t) + I < Pi x (x) q ; (t) ]2 (eqn 2.26) . 

i i 

T [ 9 R x (x)+ 0 y <p; X (x) q.(t) + 31 <Pi y (x) q.(t) ] 2 
i ' i 

+ [ X <Pi Z (x) q;(t) ] 2 
i 

= [ Z <Pi y (x) qj(t) ] 2 + [ z <Pi X (x) q ; (t) ] 2 

i i 

- 2 9 Z ( P;(x) q ; (t) Z ( Pi X (x) q ; (t) 

i i 

+ R x (x) 2 9 2 + 9 2 [ V (pj X ( x) q.(t) ] 2 + [ Z ( Pf( x ) qj(t) ] 2 

i i 

+ 2 R x (x) 9 2 q.(t) + 2 9 y <p ; x (x) q.(t) Z <Pf( x ) qj(t) 

i 

+ 2 R x (x) 9 Z ( Pi y (x) q.(t) + [ Z ( P, Z ( X ) q,(t) ] 2 
i i 

= 0 2 [ { Z <pj y (x) qj(0 ) 2 + { Z <p x ( x ) q,( t) } 2 + ( Z ( P i z (x) qj(t) } 2 ] 

i i i 

+ [ { Z <P; X (x) q ; (t) } 2 - { Z <P- y (*) q,(t) } 2 + ( Z <P; Z (x) q;(t) } 2 ] 
i i i 
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